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Role of center vortices in chiral symmetry breaking in SU(3) gauge theory 
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We study the behavior of the AsqTad quark propagator in Landau gauge on SU(3) Yang-Mills 
gauge configurations under the removal of center vortices. In SU(2) gauge theory, center vortices 
have been observed to generate chiral symmetry breaking and dominate the infrared behavior of 
the quark propagator. In contrast, we report a weak dependence on the vortex content of the gauge 
configurations, including the survival of dynamical mass generation on configurations with vanishing 
string tension. 
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I. INTRODUCTION 

The strong nuclear force has two key features: the dy- 
namical breaking of chiral symmetry (D%SB) and the 
confinement of color-charged states. It is tempting to at- 
tribute these two phenomena to a single underlying mech- 
anism, an idea supported by finite-temperature stud- 
ies where the deconfmement and chiral restoration tran- 
sitions are observed to occur at coincident tempera- 
tures Low- lying modes of the quark operator, known 
to dominate D^SB, are also correlated with the finite- 
temperature transition of the Polyakov loop, and hence 
confinement @-[Ej]. 

Over the recent past, evidence has been accumulated 
by means of lattice gauge theories that both phenomenon 
are caused by certain low energy degrees of freedoms. 
In specific gauges, these degrees of freedom appear as 
colour-magnetic monopoles [6I41] or center vortices 0- 
[Tlj . The idea that center fluxes disorder Wilson loops, 
and therefore lead to confinement, is an old one [HI, [l3[ 
and over the last couple of decades a great deal of work 
has been done in Lattice Gauge Theory on such objects, 
principally in SU(2) Yang-Mills theory. It turned out 
to be difficult to define the vortex content of Yang-Mills 
theory in a physically sensible way. It took until the late 
nineties until a successful definition was given [l4| and 
the relevance of vortices in the continuum limit was es- 
tablished The recovery of the string tension from 
"vortex-only" SU(2) gauge co nfigura tions (i.e., Z2 pro- 
jected from SU(2)) was shown [9l4ll|. the finite temper- 
ature deconfincmcnt transition was understood in terms 
of vortex properties [IrlJlq] and a connection to D%SB 
was discovered [19 - l2ll |. 

The use of Landau gauge Green's functions as probes 
of D%SB and confinement is an active area of research 
(see, e.g., [HI, [23| for a review). It is known, for exam- 
ple, gluon propagator violates spect ral p ositivity, which is 
consistent with gluon confinement [24l426l |. In the quark 



propagator the Dirac scalar part, related at large mo- 
menta to the perturbative running mass, is enhanced at 
low momenta, even in the chiral limit [27L |28| : a demon- 
stration of D^SB. One feature of this approach is that 
it allows one to make statements about light quarks, as 
opposed to the static potential of the Wilson loop. In 
SU(2) gauge theory the infrared properties of the quark 
propagator were found to be dominated by center vor- 
tices [H, H(| • Unfortunately, the vortex picture for the 
gauge group SU(3) is less clear: while vortex removal 
eliminates the linear rise of the static quark potential at 
large distances, the string tension of vortex only config- 
urations falls short by roughly a factor 2/3 [20L l3ll. |32]| . 

To gain further insights into the SU(3) vortex picture, 
we here investigate the SU(3) quark propagator under the 
removal of center vortices. We will find that mass gener- 
ation remains intact even after removing center vortices, 
while the string tension vanishes as expected. 



II. CENTER VORTICES 

We will identify center vortices in SU(3) Yang-Mills 
lattice gauge configurations using standard methods. 
Having generated gauge configurations, we will rotate 
them to Direct Maximal Center gauge then project the 
gauge links onto the nearest center element. Each config- 
uration can then be decomposed into two pieces: the cen- 
ter element and "the rest" . An appealing result of such 
a decomposition would be the identification of separate 
short and long-ranged pieces, such as seen in SU(2) gauge 
theory [29|; that is, that this decomposition corresponds 
to a separation of infrared (vortex-only) and ultraviolet 
(vortex-removed) physics. Finally, to study the propaga- 
tors, the vortex-only and vortex-removed configurations 
are rotated to Landau gauge and the quark propagators 
calculated. 

A statistical ensemble of lattice gauge configurations 
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is generated using the Liischer-Weisz [33J mean-field im- 
proved action, 



^ = ^E^ e tr(l - tf sq (x)) 



12u; 



^\lle tr(l - U Iect (x)) , 



where U sq (x) is the plaquette and U Tec t(x) denotes the 
rectangular 1x2 and 2x1 loops. For the tadpole im- 
provement factor we employ the gauge-invariant plaque- 
tte measure 



(/,] = ( -TZetT(U sq ] 



1/4 



(1) 



A. Maximal Center Gauge 

In order to identify the center fluxes of a given lattice 
configuration it is common to use gauge fixing and cen- 
ter projection. The center fluxes through an elementary 
plaquette are represented by center link elements Zp(x) 
which take values in the center group Z 3 C SU(3): 

Zp(x) = expj i^-m^x) j , m^{x) G {-1, 0, 1} . 

It is a non-trivial task to find a definition of the center 
links that is sensible in the continuum limit. The follow- 
ing definition has turned out to be fruitful [H, [III, H(| : 



(2) 



This has an intuitive interpretation: After a suitable 
gauge transformation n(x), we look for those center links 
Z^{x) that represent best a given link U fl (x). Eq. ((2|) im- 
plies that the overlap between the gauged links and the 
center links is maximized: 



5>e WtU^{x)ZI(x) 



(3) 



Hence, we will exploit the gauge degrees of freedom to 
bring U^(x) as close as possible to a center element. As- 
suming that the deviations of U^(x) from a center el- 
ement are small, one might approximately solve ([3]) by 
setting 



Z,(x) « ±TrUf!(x) , or Z^x) « [^TrU^x) 



(4) 



One gauge condition for determining the gauge transfor- 
mation is, 



max, 



(5) 



This gauge conditions specifies a particular Maximal 
Center Gauge, known in the literature as the 'mesonic' 
center gauge [|l |Hj]. 



B. Center Projection and Vortex Removal 

Once the optimal choice for the gauge transformation 
fi(x) is obtained, the center links Z^(x) are obtained from 
the gauged links U^(x) by center projection. Decompos- 
ing a particular link, 



r^(x) exp(iip M (x)) 



(6) 



where r^(x) is real and f^{x) G [— n,ir), Eq. ^ implies 
that we locally maximize 



2tt 



-m M (x) 



max 



(7) 



Hence, the integer m M (x) G { — 1,0,1} closest to 
3 ip^{x)/2n is chosen. Once the center links Z^{X) are 
obtained in this way, center fluxes ^^(x) are detected 
from the center plaquettes 



Z^x) Z v (x + /i) Zl(x + v) Zlix) 
2tt 



where 



,„(x) G {-1,0,1} 



(8) 



We say that a particular plaquette (/i, v\ x) is intersected 
by nontrivial center flux if (j)^ v (x) ^ 0. It can be shown, 
using the Z3 Bianchi identity, that the set of plaquettes 
that carry non-trivial center flux form closed surfaces on 
the dual lattice. These surfaces define the world sheets 
of Z3 vortices. The theory without center fluxes (vortex- 
removed configurations) is defined from the link elements 



U,(x) ee Uf}{x) Zlix) 



Numerical results 



(9) 



The configurations are fixed to Maximal Center Gauge 
by maximizing the gauge fixing functional ([5]) with the 
help of a local update algorithm. The algorithm is pre- 
sented in detail in (35j. Lattice sizes and simulation pa- 
rameters are listed in Table UJ 

In Fig. [TJ we show the final value of the gauge fixing 
functional i? mos for several values of the lattice spacing 
a. Generically, increasing values for i? mes are obtained 
for decreasing lattice spacing. This indicates that the 
overlap of the full configurations with pure center ones 
increases towards the continuum limit. 

If p denotes the planar vortex area density then the 
quantity pa 2 can be interpreted as the probability that a 
given plaquette carries non-trivial center flux. We have 
calculated pa 2 for several lattice spacings (see Table [TJ 
by counting the number of plaquettes with non-trivial 
center fluxes and then dividing this number by the to- 
tal number of plaquettes on the lattice. The interesting 



3 



0.76- 



0.72 



Pi 




0.68 



0.64 



1.05 0.1 0.15 0.2 0.25 

a [fin] 



0.5 



0.4 



0.3 



0.2 



0.1 



□ 12 A 3x24 lattice 
O 16 A 3x32 lattice 



0.1 0.15 0.2 0.25 

a [fin] 



FIG. 1: Left panel: The gauge fixing functional i? mes after MCG gauge fixing. Right panel: the vortex area density p in units 
of the string tension a as function of the lattice spacing a. 
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4.10 


12 3 x 24 


15 


0.611(20) 


0.272(9) 


0.1414(4) 


0.379(25) 


4.38 


16 3 x 32 


100 


0.368(5) 


0.165(3) 


0.0539(2) 


0.398(30) 


4.53 


16 3 x 32 
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0.299(11) 


0.134(5) 


0.0339(2) 


0.380(28) 


4.60 


16 3 x 32 
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0.272(11) 


0.122(5) 


0.0281(2) 


0.380(31) 


4.60 


12 3 x 24 


15 


0.272(11) 


0.122(5) 


0.0289(5) 


0.391(32) 


4.80 


16 3 x 32 


100 


0.207(5) 


0.093(2) 


0.0173(2) 


0.404(20) 



TABLE I: Simulation parameters /3, volumes, string tension 
ay/a, lattice spacings a and vortex densities. The values for 
the lattice spacings for the 16 3 x 32 lattices have been obtained 
by using 50 configurations each. For the small /3 = 4.60 lattice 
estimates are taken from the larger lattice. 



observation is that the planar vortex density, p, is inde- 
pendent of the lattice spacing a (Fig. [T]) and therefore 
has a sensible continuum interpretation. This behavior 
is in accordance with the behavior of the SU(2) vortex 
density, and confirms earlier findings for the gauge group 
SU(3) 
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FIG. 2: The quark anti-quark potential V(r) as function of 
the quark anti-quark distance r for full configurations (open 
circles), vortex-only configurations (full symbols) and vortex- 
removed configurations (open symbols). 



Now, we calculate the static quark anti-quark potential 
for ensembles with full SU(3) links, for vortex-only con- 
figurations and for vortex-removed configurations. We 
observe that the vortex-removed configurations show no 
sign of a confining potential. On the other hand, the 
vortex-only ensembles give rise to around 60% of the 
string tension. This confirms earlier findings [2fjl |3~H and 
is in sharp contrast to the case of SU(2) gauge theory: 
there, the vortices reproduce a great deal of the full string 
tension fill. 



III. QUARK PROPAGATOR ON THE LATTICE 

In a covariant gauge in the continuum, Lorentz invari- 
ance allows us to decompose the full quark propagator 
into Dirac vector and scalar pieces. In momentum space, 
the renormalized Euclidean space quark propagator has 
the form 



1 



z(C;p 2 ) 



(10) 
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where £ is the renormalization point. 

When the quark-gluon interactions are turned off, the 
quark propagator takes its tree-level form 

5(0) ^) = ^r^' ( n ) 

ip + m 

where m is the bare quark mass. When the interactions 
with the gluon field are turned on we have 

S (0)(p) _> S haie (a;p) = Z 2 ((;a)S(C; P ) , (12) 

where a is the regularization parameter (i.e., the lattice 
spacing) and ^(C; a ) is the renormalization constant. In 
the MOM scheme it is chosen so as to ensure tree-level be- 
havior at the renormalization point, Z((;( 2 ) = 1. Note 
that M(p 2 ) is renormalization point independent, i.e., 
since 5(£;p) is multiplicatively renormalizable all of the 
renormalization-point dependence is carried by Z((,;p 2 ). 
For simplicity of notation we suppress the a-dependence 
of the bare quantities. 

In this work we use the AsqTad quark action [36| be- 
cause of its excellent scaling and rotational symmetry 
properties [27l |. The Dirac scalar and vector functions, 
M(p 2 ) and Z(p 2 ) are extracted from the propagator us- 
ing the techniques described in detail in Ref. [37j . 

When analysing our results we will sometimes find it 
convenient to use a "cylinder cut" [381 ] . where we se- 
lect only data with four-momentum lying near the four- 
dimensional diagonal. This is motivated by the observa- 
tion that for a given momentum squared, (p 2 ), choosing 
the smallest momentum values of each of the Cartesian 
components, p^, should minimize finite lattice spacing ar- 
tifacts. By eliminating points most likely to be affected 
by hyper-cubic lattice artifacts it is easier to draw robust 
conclusions. 



IV. THE INFLUENCE OF CENTER VORTICES 

The Landau gauge quark propagator is calculated on 
the 16 3 x 32 configurations at = 4.60. The quark mass 
and wave- function renormalization functions of the origi- 
nal untouched gauge configurations are illustrated in Fig- 
ured Here symbols are used to identify momenta having 
a particular orientation within the lattice. Triangles de- 
note momenta lying along the Cartesian time direction 
(the long dimension), squares denote momenta oriented 
along one spatial Cartesian direction, and diamonds de- 
note momenta oriented along the lattice four diagonal. A 
comparison of triangles and squares is useful in revealing 
finite volume effects at small momenta. 

As is well-known [22|, [53] , the mass function is strongly 
enhanced in the infrared. This is true even in the chiral 
limit: a clear demonstration of dynamical chiral symme- 
try breaking. The infrared value of around 350 MeV is 
consistent with the constituent quark model. Z{q 2 ) is 
somewhat suppressed in the infrared. A study on larger 



lattices reveals a flattening of both M(q 2 ) and Z(q 2 ) be- 
low around 500 MeV This is significant for confine- 
ment, because an Euclidean propagator cannot have a 
point of inflexion and adhere to reflection positivity [28j j . 

Figure [5] shows the mass and wave- function renormal- 
ization functions after removing center vortices. Mass 
generation associated with dynamical chiral symmetry 
breaking is almost as strong after removing the center 
vortices as it was before. A roughening of the mass func- 
tion at large momenta suggests that the removal of center 
vortices introduces significant noise into the gauge field 
configurations giving rise to a larger effective mass. Z(p 2 ) 
is similarly weakly altered, being slightly noisier and hav- 
ing less infrared suppression than on the full configura- 
tions. 

Figure [5] shows a direct comparison of the quark prop- 
agator on the full and vortex-removed configurations. 
Data has been cylinder cut to facilitate a detailed com- 
parison. The wave-function renormalization function has 
been renormalized so that Z(q 2 ) — 1 at the largest mo- 
mentum considered on the lattice. Only below about 1 
GeV is there any significant difference between the full 
and vortex-removed results. It is possible that the re- 
moval of center vortices has caused Z(p 2 ) to straighten 
out, which could restore reflection positivity and hence 
be a sign of deconfinement. 

The mass function does not undergo a multiplicative 
renormalization, as described in Sec. IIIII However, re- 
moving the center vortices has significantly increased the 
running mass as displayed in the ultraviolet regime of 
the mass function. An alternative analysis is to compare 
full and vortex-removed results with bare quark masses 
adjusted to provide matched running quark masses. Fig- 
ure[B]illustratcs the persistent nature of the mass function 
under vortex removal. In this case we see that removing 
the vortices suppresses M(q 2 ) near zero four- momentum 
by about 15% compared to the full configurations, weak- 
ening - but by no means eliminating - strong infrared 
enhancement. Either way, there is still plenty of dynam- 
ical mass generation. 

To further explore the infrared nature of the quark 
mass function we turn our focus to the value of the mass 
function at the smallest nontrivial momentum available 
on our lattice, <^ in = 0.10 GeV 2 . Figure [7] compares 
the mass function at q^^ for a variety of bare quark 
masses, itiq. In the left-hand plot, Af(g I 2 nin ) is compared 
directly without an adjustment of the bare quark mass, 
whereas the right hand plot compares M(q 2 ain ) with the 
input bare masses adjusted to provide similar renormal- 
ized quark masses, m q at q = 3.0 GeV. Linear fits are 
sufficient to describe the data and indicate significant 
dynamical mass generation in the chiral limit. 

In the early days of Dyson-Schwinger studies of QCD, 
D^SB was attributed to the interaction strength in the 
quark sector provided by an effective 1-gluon exchange. 
Although it is now clear [4(| that vertex corrections play 
an important role in the quark IR sector, it is unlikely 
that a loss of gluonic interaction strength (as displayed by 
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FIG. 3: The Landau gauge quark propagator. The left panel shows the mass function M(q 2 ) and the right panel the wave- 
function renormalization function Z(q 2 ) for moa — 0.048. The infrared enhancement of the mass function demonstrates DxSB. 
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FIG. 4: Landau gauge quark propagator for moa = 0.048 following the removal of center vortices. DxSB still clearly dominates 
the mass function. Both functions are somewhat flatter than on the full configurations. 



the full gluon propagator) goes unnoticed with regard to 
DxSB. On the basis of our findings above, it is therefore 
interesting to study the effect of center vortex removal 
on the Landau gauge gluon propagator. 

As far as the full gluon propagator, D(q 2 ), is con- 
cerned, it is known to be infrared enhanced, but finite 
at zero four- momentum [4ll - |44j ]. This can be seen in 
Fig. [8] from the gluon dressing function, q 2 D(q 2 ) of the 
Landau-gauge gluon propagator. At high momenta, the 
dressing function logarithmically decreases with momen- 
tum, while it is enhanced at intermediate momenta with 
a maximum near 1 GeV. The turn over indicates a vi- 
olation of positivity, as explicitly shown in Ref. f24U26j |. 
The sa me p icture was found in full QCD with light sea 
quarks [26[ [45| . Also shown in Fig. [5] is the dressing 
function upon vortex removal. As for the gauge group 
SU(2) [46|, [47j, we find that the infrared enhancement 
is largely suppressed when center vortices are removed. 
This is particularly remarkable in the light of our previ- 



ous findings: vortex removal strongly reduces the gluonic 
interactions strength, but dynamical mass generation is 
largely unaffected. 



On a final note, we have also investigated the role 
of center vortices defined using Laplacian Center Gauge 
(LCG) [48|, where the Laplacian gauge construction re- 
moves the gauge-fixing ambiguity. While vortices defined 
this way do account for the full string tension, the vor- 
tex density diverges in the continuum limit [20(. On a 
practical side, there is an abundance of LCG vortices 
in the vacuum. Upon removing these vortices the con- 
figurations become extremely rough. The mass function 
revealed following LCG-vortex removal bares little resem- 
blance to the original mass function. It is dominated by 
noise at all distance scales as illustrated in Fig. |H1 where 
the scale has been adjusted to accommodate the results. 
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FIG. 5: Landau gauge quark propagator for m^a — 0.048. Open symbols denote the propagator obtained from the original 
gauge field configurations whereas the filled symbols denote the propagator following the removal of center vortices. Z(q 2 ) is 
renormalized to one at the largest accessible momentum point. 
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FIG. 6: The Landau gauge quark propagator with mod = 
0.048 from the original configurations (open symbols) is com- 
pared with the propagator obtained from the vortex-removed 
configurations with mod = 0.024 (filled symbols) selected to 
match the renormalized quark mass in the ultraviolet regime. 



by the persistent infrared enhancement of the quark mass 
function. The analogy to the SU(2) gauge group [2{| ends 
here: in SU(2), vortex removal implies the restoration of 
chiral symmetry [l9l - [2"lj |. 

We do stress that the key to SU(3) center vortex matter 
might not have been found yet. This is most obvious 
from the impact of the presently defined SU(3) vortices 
on the static quark potential: the string tension vanishes 
on vortex-removed configurations, but only of the order 
two thirds of the full string tension is recovered on vortex- 
only configurations. 

Whether the phenomenon of D^SB disentangles from 
quark confinement for the SU(3) gauge group (in contrast 
to the SU(2) case) or whether an improved definition of 
SU(3) vortex texture is yet to be discovered needs fur- 
ther investigation. Yet our findings offer the intriguing 
possibility to separately trace out the impact of confine- 
ment and the impact of D^SB on hadronic observables, 
at least in SU(3). 



V. DISCUSSION 

Using the SU(3) vortex picture defined by means of the 
mesonic version of the Maximal Center Gauge [H], HH 
allows us to switch to non-confining QCD upon vortex 
removal. As for the gauge group SU(2), an inspection of 
the gluonic dressing function shows a strong decline of 
gluonic interaction strength. This alone stirs the expec- 
tation that chiral symmetry might be restored as well. 

By contrast to this expectation, we find that the re- 
moval of center vortices from our configurations has done 
little to interfere with chiral symmetry breaking, as seen 
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